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We construct in the context of the AdS/CFT correspondence degenerate composite operators in 
the conformal field theory that are holographically dual to degenerate stars in anti de Sitter space. 
We calculate the effect of the gravitational back-reaction using the Tolman-Oppenheimer-Volkoff 
equations, and determine the " Chandrasekhar limit" beyond which the star undergoes gravitational 
collapse towards a black hole. 
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I. INTRODUCTION 

A lot of progress has been made in string theory in 
understanding the most extreme form of matter, namely 
black holes. String theory should also be able to de- 
scribe less extreme forms of matter, such as the degener- 
ate matter that makes up neutron stars. A neutron star 
is able to withstand the force of gravity as long as the 
total mass stays below a critical value referred to as the 
Oppenheimer-Volkoff limit pQ. This critical mass is de- 
termined by solving the Einstein equations for spherically 
symmetric pressure and energy density distributions, also 
known as the Tolman-Oppenheimer-Volkoff equations. 

ft would be especially interesting to study the collapse 
of a neutron star towards a black hole in a context with a 
holographic description. In this letter we take a first step 
in this direction by studying degenerate stars in Anti-de 
Sitter space, and constructing the corresponding states 
within the conformal field theory. One of our aims is to 
determine the OV limit for these degenerate stars and to 
give its holographic interpretation in terms of the CFT. 
Previous studies of stars in AdS are found in [2J. 

We begin in section 11 by describing degenerate com- 
posite operators in the CFT constructed out of many 
fermionic single trace primary fields. In section III we 
show that these operators correspond to a degenerate 
Fermi gas in AdS space, and we present a hydrodynamic 
description. We subsequently define a double scaling 
limit, and argue that in this limit the mass and Fermi 
energy receive corrections due to self gravity only. The 
TOV equations in AdS space and their analytic proper- 
ties are described in section IV. In section V we present 
our numerical results, compare with the boundary calcu- 
lations and conclude with remarks on the interpretation. 



II. DEGENERATE CONFORMAL OPERATORS 

We consider a CFT defined on the cylinder R x S^ 1 
with metric 



ds 1 



-dt 2 +fdfl 2 d _ 11 



(1) 



where dfl 2 { _ 1 represents the standard metric on S d 1 . 
The cylinder is identified with the boundary of global 



anti de Sitter space AdSd+i with curvature radius I. To 
make use of the operator state correspondence, we ana- 
lytically continue to euclidean time and conformally map 
the cylinder on the euclidean plane R d by 

t — > z£log |af|. 

The angular coordinates f2 are identified in the obvious 
way. States in the CFT correspond to operators via 



lim I vac) = 1$) 

1*1—0 



(2) 



We now assume that the CFT contains bosonic as well 
as fermionic fields, and has a large c limit in which one 
can distinguish single and multi trace operators [S]. 

Our aim is to find a degenerate state made out of 
fermions. So let us choose a fermionic single trace pri- 
mary operator * with conformal dimension A . * is as- 
sumed to be complex so that it differs from its conjugate 
* . We first consider the limit c — > oo, while keeping all 
other quantities fixed. In this limit single trace primaries 
behave as "free" fields, since three and higher point func- 
tions are suppressed by powers of 1/^/c, see footnote [5]. 
By spin-statistics, fermionic operators anti-commute 



(3) 



and hence they generate a free fermionic Fock space. The 
lowest dimension composites of two fields W are given by 
^d^, where di denotes the derivative with respect to x*. 
We continue with a third fermionic operator, a fourth, up 
to a very large number, each time making sure that the 
next operator is anti-symmetrized with respect to the 
previous ones by acting with one, two , three derivatives, 
etc. For a degenerate state one uses the smallest number 
of derivatives necessary to avoid that any given combina- 
tion occurs twice. For simplicity we ignore spinor indices. 
In this way we arrive at the following composite operator 



* = * ]j aw n ii didjdk* . . . 



{ii,i 2 ,— in F } 



(4) 



The fermionic primaries are arranged in "shells" labeled 
by the number of derivatives, where in each shell the 



2 



product is taken over all possible n-tuples of the deriva- 
tives di. The number of derivatives in the last shell is 
denoted by np. We consider the case with a completely 
filled last shell, so that we have a rotationally symmetric 
operator. The total number of fields contained in $ may 
then be expressed in terms of binomial coefficients as 



N 



It F 

E 

77 = 



n+d-1 
d-1 



(•5) 



Since each derivative carries conformal dimension one, 
each operator in the n-th shell has dimension n + Ao. 
Hence the total conformal dimension of <I> equals 



77=0 



A„)|"^-' 



(0) 



The primary operator \1/ corresponds to a fermionic field 
ip in the bulk. The number of derivatives acting on ^ 
is directly related to the quantum number n that de- 
termines the energy of the partial waves of ip, see e.g. 
[3]. The operator $ represents a degenerate many parti- 
cle state, whose complete wave function is given by the 
Slater determinant of all partial waves with n < np ■ 

Let us now introduce the total mass M, the fermion 
mass m, and the Fermi energy ep via 



rn 



M 



e F = {n F + A )/£. 



We are interested in the limit of large np (and large Ao) 
so it is more convenient to rewrite the sums in the ex- 
pressions §5§ and (J6j as integrals. This gives 

deg(e) , 



N(e F ) 



M(e F ) 



/ deeg(e) 



where the density of states 17(e) is given by 

t d 



ff(<0 = 



Ae-mf- 1 



(7) 



(8) 



This expression is valid in the strict large c limit. 

In the following we define a double scaling limit that 
takes np and c simultaneously to infinity, while keeping 
the ratio A/c fixed. In this limit the 't Hooft suppression 
partially breaks down, and therefore there will be correc- 
tions to the conformal dimension A and density of states 
g(e) defined in The main correction to A turns out 
to come from i-channel exchange of the stress tensor, and 
can be estimated to be of the order of A 2 /c, see [¥ . This 
corresponds to the inclusion of self gravity in the bulk. 

When operator mixing sets in, the state |<I>) has to be 
defined a little more carefully. Suppose we have a state 
\Q)n composed out of N operators 'J. We then define 



where one chooses ep so that one picks up (one of) the 
lowest order term(s) on the right hand side. For A/c — > 
this should lead to the operator Q. This procedure is 
well defined as long as the density of states g(ep) remains 
finite. When g(ep) diverges the definition of <i> becomes 
obscure. We will see this is precisely what happens at 
the Oppcnhcimer-Volkoff limit. 



III. HYDRODYNAMIC DESCRIPTION 

We return to the case where first c — > oo, and then 
take a subsequent limit in which the number of particles 
per AdS volume goes to infinity while keeping the ra- 
tio ep/m fixed. In this hydrodynamic limit the particles 
form an ideal degenerate Fermi gas. The number of par- 
ticles within one AdS radius is very large, and therefore 
the energy density p and pressure p obey the same equa- 
tion of state as in flat space. Specifically, the particle 
number density n equals the volume of the Fermi sphere 



V d . 



d(2ir) d 



Kp, 



(10) 



where V^-i denotes the area of S d 1 . The Fermi mo- 
mentum hp is related to the chemical potential p via 



k'p 



(11) 



The energy density and pressure are given by the "text- 
book" expressions 



P 



Vd-i 
(27r) d 

Vd-x 
d(2?r)' 



dk k d - 1 y / k 2 + m 2 



dk 

o 



k d+i 



and obey the standard thermodynamic relations 
dp — pdn, p + p= pn. 



(12) 



(13) 



We want to determine the density and pressure profile 
corresponding to the operator We will do so first for 
the case of pure AdS, but in the following section we gen- 
eralize these results to include self-gravity. Anticipating 
these generalizations, let us write the metric in the form 



ds 2 = -A{rfdt 2 + B{r) 2 dr 2 + ^dti 2 ^. 



(14) 



The stress energy tensor of the Fermi gas is given by the 
standard hydrodynamic expression 



T^u = {p + pyu^Uv + pg^ 



(15) 



where it M is a static velocity field: u^dx^ — A(r)dt. The 
radial profiles of p and p are determined by imposing 
stress energy conservation. This leads to the condition 



dp 1 dA 

dr + Altr {p + P) 



0. 



(16) 



This equation is surprisingly easy to solve. By making 
use of the identities (13 1 one easily verifies that (16) is 



satisfied when the chemical potential obeys 



p{r) 



A( r y 



(17) 



where at this stage ep is an arbitrary constant. We con- 
clude that the radial dependence of the chemical poten- 
tial is simply due to the gravitational redshift. We now 
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observe that the constant ep equals the chemical poten- 
tial defined with respect to the time t. This means that 
it can be identified with the Fermi energy ep of the CFT. 

When the chemical potential p is equal to the fermion 
mass m, the pressure and energy density drop to zero. 



From ( 17 1 we find that this happens at r = R for which 



A(R) = 



(18) 



This defines the edge of the star. Note that the edge is 
only present when the fermions have a finite mass. 

Up to this point our discussion holds for any metric of 
the form ( 14 1 . We now specialize to the case of pure AdS 
space, so that A 2 = 1/B 2 = 1 + r 2 /£ 2 . One then finds 
that the radius R of a free Fermi gas in AdS is determined 
by the Fermi momentum at the center: R[l = kp(Q)/m. 

As a non-trivial check on our hydrodynamic descrip- 
tion one can verify that the particle number and mass 
agree with the boundary calculations. The total particle 
number is given by the integral 



N = V d -i I drr d_1 fl(r)n(r) 



(19) 



The factor B(r) comes from the volume form on the spa- 
tial section. Similarly, one obtains the mass M by inte- 
grating the energy density. Because M is defined with 
respect to the time coordinate t one has to include a red- 
shift factor A(r). For AdS this cancels the measure factor 
B(r), and hence 



M = V d -i I drr 



d-l 



p(r). 



(20) 



In fact, this expression is equivalent to the ADM mass, 
and also holds when self gravity is included. To calculate 



N and M one has to insert the equations ( 10 1 and ( 12 1 for 



the particle and energy density and include the redshift 
effect ( 17 1 into the chemical potential. The resulting in- 



tegrals can be performed analytically and (perhaps not) 
surprisingly precisely reproduce the results ^ and 

IV. INCLUDING SELF GRAVITY 

Having established that for infinite c the operator $ 
is represented by a free degenerate Fermi gas, we now 
proceed to the double scaling limit np — > oo and c — > oo 
with fixed values for 



np 
Ao 



rn 



1. 



A 



GM 
£ d - 2 ' 



One may think this leads not only to gravitational inter- 
actions, but also possibly makes the fermions unstable. 
However, the 't Hooft suppression still seems to be suffi- 
ciently active to ensure that the Fermi gas has an infinite 
life time. The possible depletion of the particle number 
density is controlled by the Boltzman equation 



dn kp 9 
— n a, 

at m 



where a represents the total cross section associated with 
the decay. In the double scaling limit n grows like 
c d/(dH) ^ wn ji c a contains at least one factor of c _1 due to 
't Hooft suppression. So, unless a contains other factors 
of c, the life time of the fermions is estimated to be of the 
order of c 1 '^ 1 ', hence is infinite for c — > oo. The decay 
rates are possibly enhanced by energy or phase space fac- 
tors. Various checks, involving dimensional arguments, 
softness of string amplitudes and the optical theorem, 
indicate that these effects do not ruin our argument. For 
a more detailed discussion we refer to jl]. 

So let us now assume that gravity is the only force 
[§]. Except for including gravity, the fermionic matter 
will be treated exactly as in the previous section. We 
parametrize the functions A(r) and B(r) in terms of two 
new functions M(r) and x( r ) as 



A 2 (r) = e 2 *M (\ - 



C d M{r) r 2 



„d~2 



(2 



B\r) = 1- 



C d M(r) | r 2 
i 2 



„d-2 



(22) 



where Newton's constant is contained in the coefficient 

16ttG 

w 



(d-l)V d -i 

In terms of M(r) and x( r ) f ne Einstein equations read 
Af'(r) - V d - lP {r)r d -\ 

X'(r) = y d _x^(p(r)+p(r))rB 2 (r). (23) 



Together with ( 16 1 these constitute the TOV equations. 
We choose as boundary condition M(0) = 0, so that M(r) 
represents the contribution to the mass from the energy 
density inside a ball of radius r, and the total mass is 
equal to M (R). By Birkhoff 's theorem the metric outside 
the star is given by AdS-Schwarzschild. Hence, 



X (r) - 0, M(r) = M, 



for r > R. 



From ( 18 1 it follows that the radius R is determined by 

£ 2 ~~ 



1 



C d M 



R d - 



(24) 



We want to use the TOV equations to calculate the cor- 
rections to the mass, the particle number and the den- 
sity of states defined through the relations Q. In fact, 
to show that the relations ^ are still valid one needs to 
use the Einstein equations! One has 



dN 
dep 



= Vd- 



= Vd- 



dr 



.d-l 



dr ( 



drr d ~ 1 e^ 



fdM' 



Be F 

dep 

dp 

dtp 

dM 



dn dB 



dep 
dB p 



\ dei 



de, 



de F B 
_ dM 
dep 



(25) 



Here we used the identity (|13|, the relation ( |17| , the 
TOV equations ( 23 1 , and the definitions of A and B in 



(21) terms of the functions M and \ 
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FIG. 1: M as a function of log/i(0) plotted against logm for 
AdS$. The black corner is the "forbidden" region fi(0) < rn. 



M 




FIG. 2: M(ef) plotted against logm. The underlying pink 
graph is the result without gravity. The units in all figures 
are M = A/c, e F = {n F +A )/c 1 / ( - m '> and m = Aq/c 1 ^* 1 '. 



V. RESULTS AND DISCUSSION. 

The TOV equations can be integrated numerically. We 
give the results for d = 4, so for gravity in AdS$. Other 
cases are qualitatively similar. Just as in [TJ |2] we find 
that the mass M reaches a maximum as a function of the 
chemical potential at the center, see FIG 1. Beyond this 
critical value the solution becomes unphysical. Similar 
plots are found for N and tp as a function of /i(0). From 
a holographic viewpoint it is more natural to regard M 
as a function of the Fermi energy ep. This leads to the 
graph shown in FIG 2. The upper edge of the colored 



[1] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 374 
(1939), R. C. Tolman, Phys. Rev. 55, 374 (1939). 

[2] D. N. Page and K. C. Phillips, Gen. Rel. Grav. 17, 1029 
(1985), V. E. Hubeny, H. Liu and M. Rangamani, JHEP 
0701, 009 (2007), V. Vaganov, [arXiv:gr-qc/0707.0864], 
J. Hammersley, Class. Quant. Grav. 25, 205-10 (2008). 

[3] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri 



part is the OV limit, while the underlying graph is the un- 
backreacted result ^ and Q . The effect of self gravity 
on the value of M is surprisingly small. The same holds 
for N. Yet on the density of states defined as g(e) = 
dN/de the effect is significant: it diverges at the OV limit 
as shown in FIG. 3. This is a clear sign of an instability. 

Of course, one would like to derive the corrections to 
the mass from the boundary perspective. Qualitatively 
it is clear that these corrections are due to operator mix- 
ing. More importantly, we like to find the meaning of the 




FIG. 3: The density of states gCe) plotted against logm. The 
underlying graph is the result (JsJ> without back-reaction. Here 
and in FIG. 2 the black region corresponds to 6f < m. 



OV limit in the CFT. In the bulk the star becomes un- 
stable under a radial density perturbation. It would be 
interesting to interpret this mode in the CFT. Beyond the 
OV limit the star presumably collapses and forms a black 
hole. In the CFT this means that the operator $ starts to 
mix with generic states, and thermalizes. This suggests 
that the collapse is associated with a phase transition 
that turns a high density (baryonic) state into a thermal 
state (quark gluon plasma). Clearly these questions and 
the stability issue need further study. 
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